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Motivated by the heavy fermion Fermi liquid (HFFL) features observed at low-T in the pyrochlore 
LWiOa, we consider a material-specific model that includes aspects of the local quantum chemistry, 
the geometrically frustrated lattice structure, and strong correlations in a single approach. In 
particular, we show how geometrical frustration (GF) gives rise to a crossover scale, T* « J, the 
intersite (AF) exchange, below which the metallic system shows HFFL features. Our scenario is 
a specific realization of the importance of GF effects in driving HFFL behavior in L1V2O4,, and 
provides a natural understanding of various puzzling features observed experimentally. 



PACS numbers: 71.28+d,71.30+h,72.10-d 

Systems of the type AB2O4 with a magnetic B site 
consistently show strange properties. In particular, sys- 
tems with B = Ti,V,Cr are interesting because the d 
electrons occupy ti g orbitals that do not hybridize with 
the oxygen orbitals along certain directions Jl|-|5| . In the 
geometrically frustrated structure, the nearest neighbor 
exchange interaction dominates with a maximization of 
frustration. Concomitantly, only 3d bands (£29) cross the 
Fermi level (//) , leading one to have to deal with the ubiq- 
uitous strong correlations. With a divalent A-ion, the B 
site has integral valence, giving rise to Mott insulating be- 
havior with frustrated magnetism, e. g., in ZnV20±. On 
the other hand, with a monovalent A site, B is mixed- 
valent, resulting in a narrow-band metal. This combi- 
nation of frustrated magnetism and correlated metal- 
lic behavior leads to a range of complex and unusual 
manifestations, from superconductivity, for example in 
LiTiiO^, via heavy-fermion Fermi liquid (HFFL) behav- 
ior in LiViOt±, to spin-glass behavior in other systems. 

The best example in the second category above is 
LiV^Oi ||, a paramagnetic metal for T > 0.01K. The 
low-T Sommerfeld constant achieves its highest value for 
a d-electron system, 0A2J/mole/K 2 at T = 1.5K. The 
resistivity, the Woods-Saxon ratio, as well as the Wilson 
ratio all exhibit behaviors expected for prototypical rare- 
earth based HFFL metals || . Importance of geometrical 
frustration (GF) for magnetic properties is shown by the 
absence of a low-T magnetic instability, and further vin- 
dicated by neutron scattering (INS) studies M, which 
reveal a response characteristic of insulating, frustrated 
magnets. This leads one to ask: What is the role of 
GF in driving HFFL behavior? More precisely, how does 
one reconcile the HFFL behavior in thermodynamics and 
transport with a magnetic response characteristic of frus- 
trated, insulating magnets? What is the role of strong 
correlations in the 3d bands? How do these two subsys- 
tems couple to (affect) each other at low-T? 

Much attention has been devoted to these issues. LDA 
bandstructure calculations indeed show that the 3c? ti g 
bands alone cross the Fermi level M. A trigonal dis- 



tortion splits the three-fold degenerate t2 g states into the 
(lower lying) singlet A\ g with a bandwidth of leV , and an 
E g doublet with a bandwidth of 2eV in the solid. Given 
the formal d 15 state of V, the Ai g band is half-filled, 
and the E g bands are quarter-filled, leading |9) to sug- 
gestions that an effective Anderson model could be used. 
However, approaches along these lines require the rather 
ad-hoc introduction of a large intersite Kondo coupling 
to counterbalance the strong local Hund's rule coupling, 
the origin of which is unclear. Moreover, frustration ef- 
fects are not important in these pictures. On the other 
side, Fulde et al. []l0| have proposed that the V lattice 
of corner-sharing tetrahedra frustrates charge ordering 
and leads instead to isolated finite chains of 5 = 1/2 
and S — 1. The gapless fermionic spin excitations of the 
5 = 1/2 chains give the large 7 coefficient of the low-T 
specific heat. This picture needs to be extended to derive 
the HFFL properties, and to make a detailed compari- 
son with INS results. Varma |nj seeks to understand 
the HFFL properties by analyzing the energetics of the 
crossover which aS/1/2 impurity must undergo in the 
Kondo effect. Frustration does not play any role here 
either. Recently, Burdin et al. [[l2| have studied the role 
of frustration effects in driving HFFL behavior. How- 
ever, frustration effects are put in by hand, limiting a 
direct comparison with the actual material. Lacroix |I3| | 
has recently sketched the outlines of a picture for HFFL 
behavior in d-band oxides with geometrical frustration. 

Here, we study a theoretical model that is material- 
specific and includes all the relevant degrees of freedom 
and explicitly address the questions above. In particular, 
we show how the HFFL properties are reconciled with the 
magnetic response characteristic of frustrated, insulating 
magnets. The choice of the model Hamiltonian is moti- 
vated by experimental constraints 0as well as by results 
of first-principles LDA calculations |f| : 

H — H s + H b + H sb , (1) 

where H s = J^ij Jij^i-Sj is the Heisenberg-like S = 1/2 
Hamiltonian describing the localized spin degrees of free- 
dom originating from the half-filled narrow singlet A\ g 
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band. In Lil^Oi, the next-nearest neighbor coupling 
may be important, though its magnitude is not reliably 
known. 

H b = - ^ tabicl^Cjba- + h.C) + U 2j ^iaf n «4 
<ij> ia 

+ U ab ^ n iantb (2) 
i,a,b 

describes the electronic degrees of freedom for the E g 
bands. And H sb = - J H J2i, a ,b Si-(oia + &ib) describes 
the coupling of these two subsystems via a local Hund's 
rule interaction. 

In what follows, H is defined on the fully- frustrated 
pyrochlore (FFL) lattice, entailing explicit consideration 
of geometrical frustration effects. To proceed in this com- 
plicated situation, and motivated by observations, we 
make the plausible assumption: The underlying A\ g spin 
configuration does affect the charge dynamics in the E g 
band, but that there is no qualitative change in the local- 
ized magnetic response coming from the feedback effects 
due to this carrier dynamics. This assumption is justi- 
fied later. With this, we follow the following strategy: (i) 
Treat the spin correlations in the Heisenberg like model 
on the FFL within a cluster approach capable of correctly 
treating the short-range fluctuations which drive the spin 
liquid behavior as observed in INS experiments, (ii) Use 
the fact that the A\ g spins are coupled to the correlated 
electrons in the E g band(s), and modify the hopping, t ab , 
to solve the modified electronic model within the d = oo 
approximation, which is the best technique to reliably 
access dynamical effects of strong, local correlations [[l4j . 
As with the A\ g spins, the E g electrons live on the FFL, 
a fact that is indeed important for a consistent under- 
standing, as we show below. 

On the FFL, the electronic dispersion relation (or 
the spin-wave dispersion for localized spins) has four 
branches corresponding to the four-sublattice structure 
inherent in this geometry; two branches are completely 
flat over the whole Brillouin zone, while the other 
two form dispersive bands in the solid. Mathemati- 
cally, e Q (k) = 2t for a = 1,2 and, e a (k) = -2t(l ± 
y/c x c y + c y c z + c z c x ) for a = 3,4 [c a = cos(k a /2)]. No- 
tice that the free dispersion alone gives a two-band struc- 
ture in the unperturbed DOS (see Fig. |2[ with U a b = 0). 
For the Ai g spins, one replaces e(k) by J(q) and t 
by J in the above equations. From the susceptibility 
data H, a large g = 2.23 is inferred, suggesting strong 
ferromagnetic coupling between the carriers in the E g 
band and the "localized" spins in the Ai g band (large 
Jh)- So the carrier hopping rate is strongly coupled 
to, and reflects the underlying (A± g ) spin correlations: 
U(S) -> t aby /l+ < Si-Sj > /25 2 |§. 

As mentioned before, we first focus on the A\ g cor- 
relations, assuming, in accordance with LDA calcula- 
tions, that this narrow band is completely occupied by 



one localized electron per site. Following |16[|, we em- 
ploy a self-consistently embedded cluster approach that 
treats the spin correlations in one tetrahedron exactly 
and mimics the influence of the remaining tetrahedra 
by inhomogeneous, sclfconsistently determined magnetic 
fields. The resulting static spin susceptibility is shown 
in Fig. [l]. X S (T) shows a Curie- Weiss (CW)-like form at 
high-T, followed by a maximum at lower T* << J. And 
X(q, T) = g(q)x{T), consistent with 0. At very low T, 
a gap opens up in the spin excitation spectrum; by fitting 
x(q, T) to an exponential form, we extract the spin corre- 
lation length, which turns out to be weakly T-dependent 
and never exceeds a lattice spacing, in qualitative agree- 
ment with Jl?]]. The resulting physical picture is that of 
a strongly fluctuating spin system, with extreme short- 
range AF correlations and no long-range order down to 
T = 0. This agrees with the absence of magnetic order 
found experimentally, as also with the fact that the CW 
constant corresponds to S = 1/2, and that @cw < 0. 
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FIG. 1. Magnetic susceptibility of the localized spins in 
the Ai g band of the pyrochlore lattice obtained from a 
self-consistent embedded cluster approach [14] . T* marks the 
crossover temperature below which the system exhibits HFFL 
behavior. 

A potentially interesting situation now occurs: the car- 
rier hopping, as mentioned above, reflects the Ai g spin 
correlations as well as intrinsic geometric frustration ef- 
fects from the FFL structure, preventing the tendency to 
charge- and/or orbital ordering (COO). This explicitly 
realizes the mechanism for suppression of COO proposed 
by Fulde et al. fiof ; however, in our approach, it does not 
lead to formation of S = 1/2 and S = 1 rings and chains. 
The difference from the CMR manganites jll| also has a 
consistent explanation within this picture: the frustrated 
kinetic energy selfconsistcntly prevents the possibility for 
the Ai g spins to align (via "double exchange") ferromag- 
netically, in contrast to what happens in manganites |L5| . 
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Next, consider the i? ff -band electrons strongly coupled 
to the Ai g spins via a local (ferromagnetic) Hund's rule 
coupling: 



H e l = ~t ^ ( C iaa C 3ba + k.C.) + U ^ n ^a\n la \ 

+ Uab ^ n ^ n ib ~ Jh^2 S i'( Cr iQ +< J ib) ■ (3) 



where a and b refer to the doubly degenerate E g states in 
the t2 g sector, defined on the FFL. With the strong Jh, 
the "double exchange" projection transforms the problem 
to that of spinless fermions, but the hopping is modulated 
by the underlying, frustrated spin correlations in the A\ g 
band. The Hamiltonian is, 



H, 



Uj(S)(c} a Cjb + h.c) + Uab^2n ia n lb , 



(4) 



where we set U, Jh — > oo. Relabelling c a — Cj and c& = 
c±, we are left with a Hubbard- like model with a non- 
trivial hopping term: 



<ij>,a 



Cj( S )( c L c i^ + h - c ) + U *b n A n H ■ ( 5 ) 



We use the d = oo approximation to solve the 
Hubbard- like model on the FFL. This interesting prob- 
lem has non-trivial solution(s), related to the one- 
electron dispersion on the FFL. In particular, close to 
n = 1, non- Fermi liquid behavior is expected due to the 
strong scattering off the completely flat bands. Fortu- 
nately, for LiV^Oi, one deals with an almost quarter- 
filled band, rendering the flat-band singularities irrele- 
vant. The above electronic model is now solved using 
the iterated perturbation theory (IPT) at slightly less 
than quarter-filling and at finite T |Ts|| . 

Fig. U shows the evolution of the one-electron lo- 
cal spectral function, p(ui). As U a b increases, a sharp 
collective Kondo-like peak appears and gets narrower 
around /i. The appearance of upper and lower "Hub- 
bard" bands (these appear as shoulder-like features in 
the range — 2 < oj < 2 in Fig. 2), reflecting suppres- 
sion of charge fluctuations near quarter-filling, is also 
clear. Formation of heavy quasiparticles is reflected in 
the quasiparticle renormalization constant, Z(/j,), defined 

a S Z( M )= [l-^M| w=M ] _1 . 

Z([i) decreases monotonically with increasing U a b- 
That this correlated metallic state with heavy fermion 
mass (m*/m = 1/Z(fi)) is a FL is clear from the fact 
that ImS(w ~ /i) = — b(u> — /i) 2 . Importance of the 
filling is shown by the fact that charge fluctuations are 
enhanced, while spin fluctuations diminish in importance 
as n is decreased further (not shown) . Our results are in 
complete agreement with those of Imai et al. fill for a 



two-orbital Hubbard model with finite U and Jh , and the 
HFFL behavior in our model should persist when large, 
finite U and Jh are included. The computed integrated 
photoemission lincshape in the HFFL phase is shown in 
Fig. ||. Signatures of strong electronic correlations in the 
Eg band are visible as shake-up features: a broad, inco- 
herent, lower Hubbard band feature well separated from 
the narrow quasiparticle resonance. This should provide 
evidence in favor of our modelling of the E g manifold. 




FIG. 2. Local spectral density (a), and real (b) and imagi- 
nary (c) parts of the s.p. self-energy for the Hubbard model 
on the fully-frustated pyrochlore lattice for diffetent values of 
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FIG. 3. Integrated photoemission lineshapes of the Hub- 
bard model on the fully-frustated pyrochlore lattice for differ- 
ent values of U a b- 

Finally, consider the question of selfconsistency due to 
the coupling of the correlated carriers to the A\ g spin 
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fluctuations. At low T < T*, the dynamical spin sus- 
ceptibility behaves like x"(^) - (fi/T*) in the HFFL 
regime (see below). The leading-order correction to 
the self-energy of the correlated carriers is Y,(iu) n ) = 
( j h/P) E m x{i^ m )G{iu n - iVL m ). In the HFFL regime, 
G(t) ~ (1/t) for long times, and S(w) remains qualita- 
tively unchanged from its FL-like form (see above) as long 
as x"(^) — £l/T*. This analysis ceases its validity above 
T* , where the susceptibility gets increasingly dominated 
by fluctuating local moment contributions. Inserting the 
resulting into the calculation for the susceptibility 

(notice that at low T, the spin susceptibility of the lo- 
calized spins is gapped) leads again to the 0,/T* form, 
apart from prefactors. Thus, at low T < T* , selfconsis- 
tency leads only to a minor quantitative modification of 
our results. 

Let us discuss the implications of our calculation. Since 
the magnetic correlation length is T-independent below 
T* << J (corresponding to the point at which x(T) 
shows a peak, see Fig. this sets the scale below which 
the hopping, t, can be treated as constant in H&. Above 
T*, the influence of the A\ g spin correlations drive the 
system into a local-moment metallic regime, but below 
this scale, HFFL-like quasiparticles develop as shown 
above. We thus identify the low-T crossover scale seen in 
experiments with the T-dependence of the frustrated spin 
dynamics in the A\ g sector coupled to correlated E g carri- 
ers. This implies that the mechanism for HFFL behavior 
in our scenario is intimately linked to the local-moment 
magnetism of the FFL, and is drastically different from 
the conventional view, where a band of uncorrelated pcj ] 
carriers collectively screens localized moments at every 
site. Our calculation explicitly realizes the suggestions of 
Lacroix pljfl , but goes much further, showing clearly how 
short-ranged, spin-liquid-like local moment correlations 
in the A\ g sector on the FFL are related to the onset 
of HFFL behavior. In our opinion, a mean-field-like de- 
coupling cannot describe the crossover at T* adequately, 
since the spin liquid behavior is driven precisely by strong 
fluctuations beyond the mean-field picture. 

Given the HFFL metallic state below T* , the elec- 
tronic specific heat, behaves like C e i(T) ~ (T/T*), with 
a peak at T*. The local dynamical spin susceptibility 
Jx"M ~ (w/T*), so the NMR relaxation rate is Kor- 
ringa like, 1/Tj. ~ (T/JT*). The uniform spin suscepti- 
bility x — 0(1/ J) at low T. The calculation of the resis- 
tivity does not involve vertex corrections within DMFT. 
At low T, in the HFFL phase, we obtain Pdc{T) ~ 
(T/T*) 2 , implying that the Kadowaki- Woods relation is 
obeyed, as in conventional heavy fermions. The Wilson 
ratio W — T\/C e i ~ 0(T*/J), again as in the tradi- 
tional case. All these results are indeed consistent with 
thermodynamic, transport and magnetic measurements 
performed on L1V2O4. 

A detailed calculation of the dynamical spin response 
of the FFL is a fascinating problem in itself, and is inti- 



mately linked to the details of the T-dependence of var- 
ious quantities as T crosses T*. We plan to address this 
more detailed issue in a longer separate work. 

To conclude, a theoretical understanding of the vari- 
ous physical features of the frustrated, 3d, heavy fermion 
metal, LiV^O^ is proposed. In particular, we have shown 
how the HFFL thermodynamics and transport can 
be reconciled with a magnetic response characteristic of 
frustrated, insulating magnets and how the experi- 
mentally observed crossover scale (T*) is directly related 
to the T dependence of the frustrated spin dynamics in 
the Ai g sector. To our knowledge, this is the first at- 
tempt that explicitly includes specific lattice structure, 
magnetic frustration and strong correlation effects within 
a single picture. Our picture is radically different from 
conventional ones, and is a concrete realization of the 
importance of geometrical frustration effects in driving 
HFFL behavior in LiV 2 4 . 
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